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具有 未 知 参 数 的 不 确定 分 数 阶 Sprott-C 系统 的 Terminal 滑 模 同步 


闫 丽 宏 
(咸阳 师范 学 院 数 学 与 信息 科学 学 院 , 陕西 咸阳 712000) 


摘 要 : 基于 滑 模 控制 的 优良 性 能 ， 论 文 探讨 了 利用 Terminal 滑 模 控制 实现 分 数 阶 混沌 系统 的 有 限时 间 同 步 问 题 ， 给 
出 了 滑 模 控制 实现 具有 未 知 参数 和 扰动 的 分 数 阶 Sprott-C 驱动 一 响应 系统 〈 阶 数 0<a<1) 的 同步 结论 。 通 过 构造 合适 
的 滑动 模 态 曲面 ， 针 对 系统 未 知 参 数 上 界 已 知 和 未 知 两 种 情况 ， 设 计 了 合适 的 分 数 阶 控 制 器 和 参数 自 适应 率 ， 结 合 分 
数 阶 微分 方程 相关 理论 和 有 限时 间 稳 定性 定理 ， 证 明了 实现 该 系统 的 同步 控制 的 结论 ， 并 对 未 知 参 数 和 扰动 上 界 进行 
了 准确 估计 。 最 后 选取 适当 参数 ， 通 过 数值 仿真 ， 验 证 了 所 给 结论 的 有 效 性 和 可 行 性 。 
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Terminal sliding-mode synchronization of uncertain fractal Sprott-C 
chaotic system with unknown parameters 


Yan Lihong 
(School of Mathematics & Information Science, Xianyang Normal University, Xianyang Shaanxi 712000, China) 


Abstract: Basing on good performance of sliding-mode controlling, the paper investigated finite-time synchronization of fractal- 
order chaotic system by Terminal sliding-mode control strategy, and realized synchronization of drive-response chaotic Sprott- 
C System (with the fractal order a: 0<a<1) with unknown parameters and disturbance by designing property sliding-mode surface 
and fractal-order controller. Moreover, one estimates the unknown parameters and upper bound of the unknown systematic 
parameters and disturbance successfully by some appropriate adaptive update laws combining with fractal-order differential 
equations theory and finite-time stability theory. Finally, by choosing proper parameters, numerical simulations provides 
illustration to show the effectiveness and feasibility of the theoretical results . 


Key words: fractional-order chaotic system; Terminal sliding-mode control; synchronization; parameter estimation 


模 控 制 器 的 切换 项 ， 并 通过 自 适 应 算法 对 其 进行 了 自 适 应 模 糙 

调节 ， 实 现 了 非 匹 配 不 确定 性 时 变 系统 的 Terminal 滑 模 控制 。 
混沌 系统 广泛 存在 于 自然 科学 和 社会 领域 中 ， 是 一 种 呈现 庄 开 宇 等 人 0 中 针对 系统 参数 援 动 和 外 界 干扰 等 不 确定 性 医 
混乱 无 序 却 包含 内 部 规则 的 确定 性 非 线性 动力 系统 。20 多 年 前 ， 素 上 界 未 知情 况 , 实现 了 MIMO 系统 的 自 适 应 Terminal 控制 器 
Ott Grebogi 和 Yorke 成 功 将 混沌 系统 镇 定 到 不 稳定 周期 轨道 上 ， ”的 设计 。Yang 等 人 05 提 出 了 非 奇异 Terminal 滑 模 控制 方法 ， 
接着 ，Ditto 和 Roy 完成 了 混沌 控制 的 实验 ， 这 为 混沌 电路 控 。 提升 了 系统 到 达 滑 模 面 并 处 于 滑动 模 态 时 的 收敛 速度 。 


0 引言 


制 、 保 密 通 信 、 流 体 消 流 分 析 03 等 问题 的 研究 开辟 了 道路 。 近年 来 ， 分 数 阶 微 积分 得 到 了 快速 的 发 展 并 被 应 用 到 诸如 
前 , 混沌 控制 的 主要 方法 有 自 适 应 控制 9、 模糊 控 制 旬 、 最 优 。 ” 非 正常 扩散 、 信 和 号 处 理 、 量 子 力学 、 系 统 控制 等 领域 。 与 传统 
空 制 中 、 神经 网 络 控制 加、 滑 模 控 制 %9 等 智能 控制 方法 。 整数 阶 理论 相 比 ， 分 数 阶 非 线性 系统 含有 更 加 丰富 的 响应 性 能 

其 中 滑 膜 控制 近年 来 受到 了 许多 学 者 的 关注 。 在 普通 的 滑 ”和 重 棒 性 。 而 滑 模 控制 亦 具 有 响应 速度 快 、 鲁 棒 性 强 和 物理 实 
模 控 制 中 ,通常 会 选择 一 个 线性 滑 模 面 ,为 了 获得 更 好 的 性 能 ， 现 简单 等 优点 。 因 此 ， 利 用 滑 模 控制 方法 来 实现 分 数 阶 混沌 系 
Xu 等 人 提出 了 Terminal 滑 模 控制 策略 中 ， 可 使 系统 在 有 限时 ” 统 的 同步 , 将 极 大 地 提高 控制 效率 。 基于 此 , 文献 [16] 研 究 了 含 


nk 


间 内 实现 对 期 望 状态 的 完全 跟踪 。 文 献 [12] 探 讨 了 非 奇异 ” 有 建 模 不 确定 参数 和 外 部 扰动 的 整数 阶 混沌 系统 的 滑 膜 控制 问 
Terminal 滑 模 控 制 器 的 设计 问题 ， 实 现 了 含有 不 确定 项 的 混沌 。 题 。 邓 立 为 等 人 [17] 利 用 输出 反馈 滑 模 控 制 实现 了 分 数 阶 超 混 
系统 的 渐 近 同步 。Tao 等 人 [3 采用 模糊 规则 设计 了 Terminal 滑 。 沌 系统 的 同步 。 文 献 [4,5,8,19] 对 分 数 阶 ( 超 ) 混 沌 系统 的 自 适 应 
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录用 稿 


滑 模 控制 问题 进行 了 分 析 研 究 。 李 特等 人 RV 设计 了 一 种 滑 模 控 
制 器 , 利用 分 数 阶 微 积分 理论 ， 实现 了 两 个 不 确定 分 数 阶 混沌 
系统 的 同步 控制 ， 并 辨识 出 了 误差 系统 中 不 确定 项 及 外 界 干 扰 
项 的 边界 ,文献 [21] 基于 lyapunov 稳定 性 理论 和 分 数 阶 微 积分 
的 Laplace 变换 , 研究 了 分 数 阶 Genesio-Tesi 系统 的 滑 模 混沌 同 
步 。 高 俊 山 等 人 [3 结合 自 适 应 控制 理论 和 滑 模 控制 理论 ， 提出 
较 传 统 滑 模 面具 有 更 快 收敛 速度 的 Terminal 滑 模 面 ,设计 的 自 
适应 滑 模 控制 律 保证 了 系统 在 有 限时 间 内 收敛 到 滑 模 面 。 由 此 
可 见 ， 利 用 滑 模 控制 实现 分 数 阶 混沌 系统 的 同步 具有 重要 的 研 
究 价 值 。 

对 于 实际 系统 ， 由 于 建 模 误差 和 工作 环境 等 因素 的 影响 ， 
随机 和 外 部 扰动 是 客观 存在 的 。 综 合 以 上 考虑 ， 本 文 将 提出 一 
种 新 的 分 数 阶 滑 模 面 ， 并 设计 Terminal 滑 模 同 步 算 法 ， 证 明 误 
差 动 态 系统 可 在 有 限时 间 内 收敛 到 平衡 点 ， 从 而 实现 含有 未 知 
参数 和 扰动 的 分 数 阶 Sprott-C 驱动 -响应 混沌 系统 的 同步 , 同时 
在 线 估 计 未 知 参 数 和 扰动 上 界 。 最 后 利用 MATLAB 数值 仿真 
验证 所 给 同步 方案 的 正确 性 和 有 效 性 。 


1 ”模型 建立 和 知识 预备 


CP3 混 沌 驱动 系统 : 


Dr” x =a(x,—X) 
Dr xX, = —a3%, 一 为 罗 (1) 
一 
其 中 : 为 ,2 ,2 分 别 为 该 系统 的 状态 变量 , 分 数 阶 次 C 满足 : 
0<C<1。 

选取 含有 未 知 
统 为 


界 随 机 扰动 和 外 部 扰动 的 分 数 阶 响应 系 


Dr y= — ytH(Y) + Ot) + ult) 
有 Dr 六 =-0 六 -ptpO)+wOD+OOD  O) 
Dr” y, = 六 —A,+n(y) + 0,(t) +u,(7) 


其 中 ，4 是 未 知 参数 a 的 估计 值 。y = (yy, 3), 1(y)， 
(7), 1 = 二 1,2,3 是 满足 假设 1 的 具有 未 知 上 界 的 随机 扰动 和 


外 部 扰动 项 ，u,(?) 为 滑 模 同步 控制 器 。 对 于 不 确定 项 和 外 部 
扰动 ， 有 如 下 假设 : 
假设 1 “假设 存在 未 知 有 界 正 数 广 ，C, ， 使 得 随机 扰动 和 
外 部 扰动 项 分 别 满足 : 
[ws pa 
leo so, 


是 合理 的 ， 


其 中 : 表示 欧 氏 范 数 。 这 一 假设 在 实践 应 用 
为 实际 发 生 的 信号 均 是 有 界 的 。 
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假设 六 中 分 别 是 未 知 随机 和 外 部 扰动 上 界 7,d; 的 估计 值 ， 
且 满 足 如 下 的 自 适 应 率 : 


Si|， (0) 一 fo, 


六 = 大 
& =k,|s,|, %(0) = 四， 


13 (4) 


记 玉 = 一 4,@ = 心 一 ,i=1,2,3 分 别 为 上 界 ,0 的 
估计 误差 。 

驱动 一 响应 系统 式 (1) 和 (2) 之 间 的 误差 为 

e(1) = y(1) —x(?) (5) 

则 系统 式 (1) 和 (2) 之 间 的 同步 问题 等 价 于 误差 系统 式 (5) 的 稳定 
性 问题 。 当 对 响应 系统 施加 合适 的 滑 模 控制 器 后 ， 误 差 e(1) 在 
零点 稳定 即 意 味 着 系统 式 () 和 (2) 之 间 实 现 滑 模 同步 。 故 可 将 同 
步 问题 转换 为 稳定 性 问题 来 研究 。 

对 误差 动态 方程 (5) 求 导 可 得 

Dr el(t) = Dr y(t)— Dr x(?) (6) 


即 


Dre =A(e,—e)+a (x —xX)+tn(y) 
+ (1) + u(t) 
Dre, = -bse — xX, — yy3tX Xs tr (y) 0) 
+ 0,(f)+u, (1) 
Dre,=(y,+X,)e,—d,+n(y)+0,(7) 
+ us(1) 


为 了 讨论 方便 ， 引 入 如 下 假设 和 引 理 : 
假设 2 假设 未 知 参数 a 满足 约束 条 件 : 


[ola 1=12.3 


其 中 区 是 正常 数 。 
引 理 1 有 限时 间 稳定 性 定理 2q。 假 设 存在 连续 正定 的 函 
数 V(1) 满足 下 面 的 微分 不 等 式 
V(D) <-eV "(1),vt>t,V()>0, 
其 中 : E> 0,0< w <1, 那么 对 于 任意 给 定 的 th，V(7) 满足 下 
面 的 微分 不 等 式 : 
VD EV HN) -elt), to stst,. 妥 
V1 之 fi,V(?) 三 0 其 中 有 限时 间 1i 为 


V (0) 
”0 ed-o0 


2 分数 阶 Sprott-C 系统 的 Terminal 滑 模 控 制 


2.1 有 限时 间 Terminal 滑 模 控制 

参照 文献 [19] 中 对 传统 整数 阶 滑 模 面 的 设计 ， 本 文 进行 了 
改进 ， 将 其 推广 到 分 数 阶 情形 ， 设 计 如 下 形式 的 Terminal 滑 
模 面 : 


录用 


其 中 : 0<ag<1, 0<0<1, 7 > 0 为 增 


稿 


5,(1) = DY le (1) + 


0 (8) 
中 sgn(D e 


多 eo dv 


益 系 数 ， 记 7 = min{ 


11772;773} 。 对 该 式 求 导 可 得 
$(D=Dre (D+ sgn(D ep ed) ©) 
当 误 差 系统 在 滑 模 面 运动 时 ， 滑 模 面 及 其 导数 需要 满足 可 


达 怕 


在 有 限时 间 i; 内 达到 


E 条 件 S (D =0, s,(7)=0。 
定理 1 在 Terminal 滑 模 面 式 (6) 上 ,误差 系统 式 (5) 的 轨迹 
平衡 点 ， 其 中 : 


Th, < io ‘(Foe o 门 J-o “71 (10) 


i=] 


一 证 明 误差 系统 满足 滑 模 面 方程 : 5,(1)=0,5,(1)=0。 
于 是 可 得 如 下 导数 : 
Dre() =-n,sgn(D" ep" e 0)| (11) 
取 Lyapunov 函数 为 
1 3 
VCD) = 2 Derle 的 
1i=1 
沿 着 时 间 f 求 导 可 得 
. 3 
VD= >, De(n) Dre,(t) = 
的 | 
1 1 1 0 
> De .CC7sgn(D7 et) Dp en)| )= 
i 
2 | 1+0 
-2n|D; RO) 
i=] 
利用 不 等 关系 式 (| ++, 站 < 
+ 
1+0 
, 1+0 1 生 7 
VD)<-27 "人 (oreo | . (12) 
i=1 


tm 内 达到 平衡 点 ， 


2.2 


节 将 讨论 误差 系统 在 滑 模 


根据 有 限时 间 稳 定性 定理 ， 可 知 误差 系统 将 在 有 限时 间 
1 满足 以 下 条 件 : 


fp < 六 :So ‘6,(10)) 小 ， or 


这 意味 着 误差 系统 (7) 在 有 限时 间 内 到 达 滑 模 面 。 
模 控制 
保证 了 误差 系统 可 在 有 限时 间 内 达到 平衡 点 。 本 
而 上 的 稳定 性 问题 。 将 驱动 系统 式 (]) 


Terminal 滑 
上 述 定 至 
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| V 瑟 人 FE 基本 
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的 未 知 参 数 上 界 ,i =1,2,3 分 为 已 知 和 未 知 两 种 情况 。 针 对 不 


同情 形 ， 将 分 别 设计 合适 的 控制 器 ， 实 现 分 数 阶 驱动 -响应 
Sprott-C 系统 的 滑 模 同步 控制 。 

情况 1 了 驱动 系统 的 参数 上 界 避 ,i 二 1,2,3 已 知 

设计 未 知 参数 4 的 自 适应 更 新 率 为 ; 


a 硅谷 (a 一 六 )8 一 和 sh)， 0(0) = Ao; 
(13) 


0 =c, (ss -6,|s,|), a,(0) = do; 


a, =c, (Xs -人 |s;|), 4,(0) = 0 


其 中 参数 c;，i = 1, 2,3 均 是 大 于 零 的 常数 ， 可 以 调节 同步 
的 收敛 速度 。 并 且 估计 误差 分 别 记 为 : 6 = 一 ,i=1, 2,3。 
定理 2 假设 驱动 系统 未 知 参数 的 上 界 蔬 ,i=1,2,3 已 知 ， 则 
在 未 知 参数 和 扰动 自 适应 率 式 (4)(13) 和 控制 率 式 (14) 作 用 下 ， 
驱动 -响应 系统 式 (0 和 (2) 可 实现 滑 模 同步 控制 。 控 制 率 如 下 


入 CQ 一 QO= 9 
u(t) = -A(e, -ei) -mn sgn( Dr e)|D, er(b| 


一 人 + +a? )sgn(s,) 


二 = C 一 
oOD = bes + yy -sen( Dr elD eo))| (14) 


一 人 + 外 + 五)sgn(s )— Xx 


万 和 1 2 
ex(D| 


友人 人 = 一 (六 十 各)e —7sgn( DY 2) 


-人 十 心 十 而 ?)sgn(s) 


构造 Lyapunov 能 量 函 数 如 下 : 


并 沿 着 时 间 f 求 导 可 得 
VOD=》 ss |[t Ft 
小 


将 参数 更 新 率 式 (13)、 扰 动 


+15 i 
0, i 


适应 率 式 (4) 分 别 代入 上 起 可 


得 


Dere oh 


7O=> (Dre,(t)+7, sgn(D™'e,(t)) 


Dy 一 .大 ;|i | + (Xs, 一 人 js,|) 
+ (a —%b)s 一 i —6, [sl) 
三 5 GAC —X)+n(y) +0()) 


+ -OO Done of sgn(D”'e, GO 


—e)+A (% 


@ tN — Vy tAX tn(y)+ 0@,(7)) 


De, of sgn(D"e,)] 


2 (6 
+5,: (worm 
二 53 ((> 十 六 )e， -6 +n(y)+0,(t) +us(t)) 
D, [sl)+ G(x)s -alsl) 


5,|)+ t, (xs, 一 0 js,|) 


Sy 


录用 稿 


0 
Dr 'e,(D)| sgn(D7 es) 


十 53 773 
将 控制 器 表达 式 (14) 代 入 上 式 有 
VOD=50 (0)+oD-(i 


+5,(D(n(y) + oD (S++ )sen(s,)) 


+ +a’ )sgn(s,)) 


+53(0) (50) + oD) (B+ +E )sen(s,)) 


引 | 
i=] i 
3 本 
> 2 256 中 
i=1 i=1 
于 和 -0=0.00 -a,)=246, -A <246,， 即 


VD < s a |s, 
<- |. 


i=l 


于 >0， 所 以 V(t) <0。 根据 LaSalle 不 变 原理 可 知 当 
1 一 +oo 时 ，si(GD 下 0。 进一步 , 根据 分 数 阶 微 积分 理论 可 知 
系统 误差 e;(1) 一 0,7 =12,3, 这 意味 着 驱动 响应 系统 式 (1) 和 (2) 
实现 了 Terminal 滑 模 同步 。 

情况 2 驱动 系统 的 参数 上 界 忆 ,i 二 1,2,3 未 知 

此 时 未 知 参数 a 的 自 适 应 更 新 率 可 设计 为 


a =c (XO—%)s, A(0) = 0; 
bh = C053, O(0) = 6; (15) 


0 = C3%5,, (0) = Oo. 


参数 c; 与 自 适 应 率 (13) 中 的 含义 一 致 。 


定理 3 在 未 知 参数 和 扰动 自 适应 率 式 (4)(15) 和 控制 率 式 
(16) 作 用 下 ， 了 驱动 一 响应 系统 式 (1) 和 (2) 可 实现 滑 模 同步 控制 。 
其 中 控制 率 如 下 : 


各 QQ 一] 一 9 
w(D) =-A (es —e) -nsgn( De) Dr re 


一 人 + +l)sgn(s,) 


入 Q— Q— 0 
W(t)= 6, + yy -1 sgn(D, ‘6,) D, ‘(0)| 


(10) 


—(KB 1+, +h )sgn(s,) — Xx 


0 
pre,()| 


ID =—(y, +%)e, 一 7 Sgn(D2” e) 


—(B+O, +h )sgn(s,) 
其 中 增益 系数 1 > 0，i=1,2,3。 

证 明 构造 和 定理 2 相同 的 Lyapunov 能 量 函数 V(1)， 
沿 时 间 f 求 导 。 然 后 分 别 将 扰动 自 适应 率 式 (4)、 参 数 更 新 率 式 
(15)、 滑 模 面 方 程 导数 式 (9) 和 控制 器 式 (16) 代 入 V(7) 可 得 
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VO=50 (nO) + oD (i+ + )sen(s)) 
+5,(D) (nO) + 0D) (b+, +1 )sgn(s,)) 
ts5(D (HO) + - (B+ +L )sgnGss)) 

lsl+ ls,) — or 


1=1 


< -> Si Sgn(s,) 
i=1 


令 1=min{1,b,4} 有 V(?)< Dssen(s)=— >ils. 不 


难 发 现 ，V(1)<0， 即 当 t > 0 时 ， Lie 数 单调 递减 。 
在 区 间 [0,t] 上 关于 时 间 t 进行 积分 ， 结 合 比 较 原理 可 得 


[ > (OE < -了 VO)de = (VO) -VD) 


< (VO -VC)) < TY < +oo 


结合 Barbalat 引 理 , 可 知 e(D) 一 0, [=12,3。 这 意味 着 驱 
动 系统 未 知 参数 上 界 未 知 时 可 实现 系统 式 ) 和 (2) 的 滑 模 同步 
控制 。 


3 ”数值 仿真 


本 章 中 选取 合适 的 参数 对 上 述 讨论 结果 进行 数值 模拟 。 为 
避免 滑 模 控 制 中 的 拌 振 现象 ， 在 实际 工程 应 用 中 ， 通 常 有 两 种 
方法 四 可 以 消除 抖 振 现象 。 第 一 种 是 将 控制 器 分 为 连续 和 切换 
两 种 ， 消 除 切换 振幅 ， 第 二 种 是 利用 高 增益 方法 ， 用 饱和 函数 
等 函数 代替 Sgn 项 。 在 仿真 实验 中 ， 本 文采 用 第 二 种 方法 ， 
参照 文献 []， 利 用 s/(|s|+e) (其 中 E=0.001) 代 蔡 
sgn(s,(1)) 项 ， 通 过 仿真 发 现 ， 可 达到 良好 的 控制 效果 。 

当 分 数 阶 次 w=0.97,a =10,m =100,4,=-1.7 时 ， 
Sprott-C 系统 是 混沌 的 。 且 有 两 个 平衡 点 Ei=(10, 10,1.7) 和 E2=(- 
10,-10, 1.7)。 混 沌 吸引 子 如 图 1 所 示 。 为 了 便于 模拟 ， 取 系统 
中 的 不 确定 项 分 别 为 nn(y)=0.2sin(37)y,,n(y) 
二 0.15sin(27)y,n(y)=0.3cos(yy) ， 外 部 扰动 是 @(7) 
二 0.2sint,@y(1)=0.3cost, C3(1)=0.4sint 。 仿 真 中 ， 扰 动 初始 
值 分 别 取 为 A0) = (0.1,0.2,0.3),Q(0) = (0.7,0.5, 0.3) 。 不 失 
一 般 性 ， 取 扰动 上 界 的 调节 参数 =1.1 态 =0.9 。 滑 模 面 (8) 中 
参数 分 别 为 (07,7p,7])=(0.5,0.5,0.4),0= 0.2 。 

当 驱 动 系统 的 参数 上 界 志 ,i =1,2,3 已 知 时 ， 上 界 取 为 
吉 =10,a,=100,a,=1.7 。 系 统 初 值 为 : x = {0.8, 一 0.2, 一 
1.9), yo = (1.6,1-4) 。 未 知 参 数 的 初 值 分 别 选 为 : (Qi 

0 ,0.) = (7.41,97.9, 一 1.92) ， 参 数 自 适 应 率 系数 选取 
为 ，c;=1, i=1,2,3 。 利 用 预 估 -校正 数值 算法 ， 对 系统 的 同步 


真 ， 结 果 如 图 1~3 所 示 。 


进行 仿 
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图 1 分 数 阶 Sprott-C 混沌 吸引 子 


time(s) 


图 2 了 驱动 响应 系统 (1) 和 (2) 的 误差 演化 (情况 1) 


0 1 2 3 4 5 6 7 


to 
© 已 
% So 


0 1 党 3 4 5 6 了 
time(s) 
20 
中 0 
-20 
0 1 3 5 6 


time(s) 


图 3 驱动 系统 (1) 的 参数 估计 曲线 @ =0.97 (情况 1) 

图 1 是 分 数 阶 w = 0.97 时 Sprott-C 系统 的 混沌 吸引 子 
当 对 驱动 响应 系统 式 (1) 和 (2) 施 加 控制 器 式 (14) 时 ， 系 统 的 误差 
如 图 2 所 示 。 从 图 2 可 以 看 到 , 系统 误差 在 有 限时 间 内 趋 于 零 ， 
这 意味 着 驱动 一 响应 网 络 式 (1) 和 (2) 实 现 了 Terminal 滑 膜 同步 。 
同时 ， 驱 动 系统 的 未 知 参数 也 得 到 了 准确 估计 ， 结 果 如 图 3 所 
示 。 可 以 看 到 ， 未 知 参数 均 收敛 到 其 精确 值 。 

事实 上 ， 当 分 数 阶 次 =0.98 时 ，Sprott-C 系统 仍然 是 混 
沌 的 ， 此 时 系统 的 混沌 吸引 子 与 两 个 不 稳定 平衡 点 B1 和 Bs 共 
存 。 选 定 初 值 = (0.5,0.3, 一 1.2), y = (3,1.5,-5 ) ， 驱 动 系 
统 未 知 参数 的 初 值 分 别 选 为 (2 ,0 ,0,) = 
(10.57,99.7,-0.53) ， 其 余 参 数 取 值 均 和 上 述 情况 相同 。 
进行 仿真 ， 此 时 驱动 系统 误差 演化 结果 如 图 4 所 示 ， 同 时 系统 
的 三 个 未 知 参数 的 估计 结果 如 图 5 所 示 。 


丽 
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time(s) 


图 4 ”驱动 响应 系统 (1) 和 (2) 的 误差 演化 (情况 1) 


< 10 站 
5 
0 2 4 6 8 10 
time(s) 
100.5 
本 De 
99.5 
0 2 4 6 8 10 
time(s) 
2 
0 
| 
Sd 
-6 
-8 
0 2 4 6 8 10 


图 5 驱动 系统 (1) 的 参数 估计 曲线 w=0.98 (情况 1) 

对 于 情况 2， 当 未 知 参数 的 上 界 未 知 时 ， 利 用 控制 器 (16) 
可 实现 系统 (1) 和 (2) 的 同步 。 固 定 取 如 情况 1 的 相关 参数 ， 反 
馈 参数 为 ，11=0.2, 4,=0.15, /=0.2 ， 初 值 取 为 2 = 
(2.6,1.2, 一 3), yo = (1 0.2,0.2) ， 未 知 参数 的 初 值 分 别 选 
为 : (6 0. 6) =(12.3,102.8,—3.16),c, =0.1,i= 
1,2,3 。 系 统 阶 次 w= 0.97 ， 误 差 演 化 曲线 和 未 知 参数 估计 分 
别 如 图 6 和 7 所 示 。 


图 6 系统 趋 于 平衡 点 El (情况 2) 


time(s) 
图 7 驱动 系统 (1) 的 参数 估计 曲线 (情况 2) 
8 是 系统 逐渐 趋 于 平衡 点 Ei 的 演化 曲线 。 同时, 经 过 计 
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算 ， 可 知 上 述 两 种 情况 中 ， 仿 真 中 实现 同步 的 时 间 与 定理 1 中 
的 理论 计算 时 间 (10) 相 吻合 。 
20r . 
< 导 10 
0 
-10 
0 2 4 6 8 
time(s) 
103 
102 
‘S 101 
100 
0 2 4 6 8 
timel(s) 
20 
8 0 | 
-20 ] 
0 2 4 6 8 
time(s) 


对 


8 ”驱动 响应 系统 (1) 和 (2) 的 误差 演化 (情况 2) 


对 


9 系统 趋 于 平衡 点 E2 (情况 2) 

固定 反馈 参数 和 系统 初 值 等 相关 参数 ， 未 知 参数 的 初 值 分 
别 选 为 (6 6)=(4.56,95.6,-2.15),c =2,i 
= 1,2,3 。 系统 可 以 稳定 到 平衡 点 E2， 仿 真 结果 可 见 图 9。 此 
时 系统 未 知 参数 的 估计 曲线 和 误差 演化 曲线 如 图 10 和 11 所 示 。 


二 we 一 
-100 
0 2 4 6 8 
time(s) 
2 > ne 
全 
95 
0 2 4 6 8 
timel(s) 
60 
40 上 
二 2 
-20 
-40 


0 刘 4 6 8 
time(s) 


图 10 驱动 系统 (1) 的 参数 估计 


线 (情况 2) 


time(s) 


对 


11 驱动 响应 系统 (1) 和 (2) 的 误差 演化 (情况 2) 
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通过 以 上 数值 仿真 可 以 看 到 ， 未 知 参数 和 随机 (外 部 ) 扰 动 
的 估计 值 均 趋 向 其 真实 值 , 驱动 系统 (1) 和 响应 系统 (2) 达 到 了 有 
限时 间 同 步 ， 充 分 说 明了 所 设计 的 同步 控制 器 和 未 知 参数 辨识 
规则 是 有 效 的 。 


4 ”结束 语 


tr 


研究 混沌 控制 的 目的 是 希望 在 更 大 范围 内 更 好 地 操纵 非 线 
性 系统 的 动态 ， 以 便 为 人 们 所 用 。 本 文 利用 分 数 阶 控制 器 对 分 
数 阶 被 控 对 象 进行 了 同步 控制 ， 探 讨 了 含有 未 知 参数 和 有 界 扰 
动 的 不 确定 分 数 阶 Sprott-C 系统 的 Terminal 滑 模 控 制 问题 ， 设 
计 了 滑 模 动态 面 函数 , 实现 了 所 构造 驱动 -响应 系统 的 滑 模 同步 
控制 ， 同 时 对 未 知 参 数 进 行 了 估计 ， 最 后 通过 数值 仿真 验证 了 
结论 的 有 效 性 。 本 文 的 方法 和 结论 具有 一 定 普 适 性 ， 可 为 同类 
型 分 数 阶 混沌 系统 的 同步 控制 研究 提供 一 定 的 参考 。 
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